Generation of Bright Soliton through the Interaction of Black Solitons 
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We report on the possibility of having two black solitons interacting inside a silica fiber that presents 
normal group-velocity dispersion, to generate a pair of solitons, a vector soliton of the black-bright 
type. The model obeys a pair of coupled nonlinear Schrodinger equations, that follows in accordance 
with a Ginzburg-Landau equation describing the anisotropic XY model. We solve the coupled equa- 
tions using a trial-orbit method, which plays a significant role when the Schrodinger equations are 
reduced to first order differential equations. 
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Solitary waves and solitons are relevant to a variety 
of nonlinear physical processes. They have made per- 
haps the geatest experimental impact in the field of non- 
linear fiber optics, which are widely known to support 
solitons of several different types |I|-q]> For instance, in 
silica fibers there are optical solitons of the bright and 
dark types. They appear because the Kerr nonlincar- 
ity in silica is always described by a positive coefficient, 
so that the sign of the group-velocity dispersion distin- 
guishes two different types of solitons in silica fibers: the 
bright soliton, for negative sign of the group- velocity dis- 
persion, and the dark soliton, for positive group-velocity 
dispersion. These two types of optical solitons were first 
reported in Ref. fl, in the case of bright solitons, and in 
Ref. @, in the case of dark solitons. They are solutions 
of the cubic nonlinear Schrodinger equation, and they 
are characterized by having distinct profile: the bright 
soliton has vanishing asymptotic behavior, so they en- 
gender trivial topological behavior. On the other hand, 
dark solitons are characterized by non vanishing asymp- 
totic behavior, so the presence of topological profile is 
one important feature of dark solitons. The central char- 
acteristic of dark solitons is the presence of a dip in its 
center, which distinguishes two possibilities: black soli- 
tons, in the case the dip reaches the botton, that is, if 
the dip goes completely to zero, and gray solitons, in the 
other cases. 

Although both bright and dark solitons are widely 
known to be present in fibers, they are not exclusive of 
fibers. They can also appear in other systems, in particu- 
lar in Bose-Einstein condensates (BEC), and in magnetic 
materials. The old prediction that an ideal gas of identi- 
cal bosons may condensate into a macroscopic quantum 
state has recently been accomplished: BECs of differ- 
ent atomic elements were first produced in ^-j^] , and are 
reviewed for instance in Ref. || . Soom after this achiev- 
ment, the presence of bright solitons in a BEC was first 
reported in Ref. Jl(j] , in the study of the dynamics of ul- 
tracold bosonic atoms loaded in an optical lattice, that 
is, in a lattice induced by the interference of an array 
of laser beams pi. The presence of dark solitons was 



reported more recently in Refs. 12 13|]. The presence 
of bright or dark solitons in BEC is directly related to 
the physics of attractive or repulsive atomic interaction, 
respectivelly, in a mean-field description that follows in 
accordance with the Gross-Pitaevskii equation |T^] . In 
magnetic films one finds dark solitons jl5|-|l7j in the form 
of a microwave magnetic envelope (MMEjT A good ex- 
ample of this was recently reported in Ref. JL8[, where 
one can find an interesting way of generating dark MME 
solitons, which opens a new route for the investigation of 
dark solitons in magnetic systems. 

In silica fibers both black and bright solitons spring as 
solutions of the cubic non-linear Schrodinger equation, 
with just a change in the sign of the term that controls the 
group-velocity dispersion in that equation. We illustrate 
the two possibilities considering the normalized equation 
that describes the electric field envelope U = U(x,z) 
inside the optical fiber 



dU± ± ld 2 U± 



\u±\ 2 u ± = o 



(1) 



dz 2 dx 2 

The variables x and z describe the reduced time and 
space coordinates inside the fiber, respectively. This 
equation takes into accound the temporal or group- 
velocity dispersion, the second term in the above 
Schrodinger equation, and weak Kerr nonlinearity, the 
third term in Eq. ([l]). The two signs represent the two 
standard situations, that generate the black and bright 
solitons. 

The presence of solitons follows after writing 
U±(x,t) = u±(x)exp(iz). We get 



ld 2 u± 



— ± u± T «j 



(2) 

(3) 
(4) 

where x is an arbitrary point, standing for the center of 
the soliton. Here u+(x) and U-{x) represent the fun- 
damental solitons, the bright and dark solitons, respec- 
tively. In fact, u-{x) is a black soliton since the dip in 



2 dx 2 

These equations are solved by 

u + (x) = ±sech(x — x) 
u_ (x) — ± tanh(x — x) 
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u 2 _{x) goes to zero, reaching the botton of the electric 
field envelope. 

We consider the possibility of describing a more gen- 
eral situation. We think of using a single fiber, an 
optical medium with normal group-velocity dispersion. 
However, we illuminate the fiber with two distinct laser 
beams, having different amplitudes and opposite circu- 
lar polarizations. The electric fields are characterized by 
the envelopes U(x, z) and V(x, z), and the system is de- 
scribed by the two normalized equations 



.dU _ l d 2 U 
dz 

;— 

dz 



F(I)U = 



2 dx 2 
1 d 2 V 

-^ + G(I)V = 



(5) 
(G) 



We use / = \U\ 2 + \V\ 2 to describe the total intensity 
inside the fiber, and now the envelopes of two interacting 
laser beams inside the optical medium are characterized 
by the functions F(I) and G(I), which respond to non- 
linear Kerr interactions inside the silica fiber. 

We describe the case of two different laser beams, that 
interact inside the silica fiber with the envelope of the 
electric fields in the form 



U(x, z) — u(x) e l 
V(x,z) = v(x) e k 



(7) 

(8) 



Here r is real and positive, and parametrizes the rela- 
tive propagation constant of the vector constituents. We 
substitute Eqs. (]?]) and (^) into Eqs.(||) and (||) to obtain 



1 d 2 u 
2dx 2 
lfy 
2dx~ 



— u + F r (u, v) u 



o3Z2 = -rv + G r (u,v) v 



(9) 
(10) 



We normalize the system in a way such that, in the ab- 
sence of v we get F r (u, 0) = u 2 , and in the absence of u 
we get G r (0, v) — r 2 v 2 . In the first case, for v — we get 



1 d 2 u 
^dx 2 

whose solutions are 



-u + u° 



Ub{x) = ± tanh (x — x) 
The second case is similar, leading to 



v b (x) — tanhfVr ( 



(11) 



(12) 



(13) 



These solutions represent black soliton solutions. The 
system we are preparing is such that in the reduced case 
of just one component, obtained after turning off one of 
the two laser beams, the interactions between the optical 
fiber and the beam give rise to a black soliton. Each one 
of the black solitons has its own characteristics, which 



are naturally related to each one of the two laser beams 
since we are not changing the optical medium, the fiber. 
The individual black solitons are different, giving rise to 
an asymmetrical arrangement that will allow interesting 
novelties. The idea is similar to those in Refs. [p~9[~pl| , 
and the results remind us of other effects: in quantum op- 
tics, a superposition of two (antisqueezed) number states 
may exhibit squeezing effect [ p2| , a property not shown 
by any isolated component; also, a superposition of two 
(most classical) coherent states may result in a nonclassi- 
cal (Schrodinger's cat) state, as observed experimentally 
for light fields [^3) and for trapped ions p4| . 

The main characteristics of the two black solitons are 
depicted in Fig. 1, where we plot both the Ub(x) and 
Vb(x). We notice that the amplitude (a) and width (w) 
of Ub(x) are such that a u = w u = 1, but for Vb(x) they 
become a v = w v = \J\jr. Thus, in the following we will 
consider the case re (0, 1). 
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FIG. 1. The two black solitons, ul(x) and v%(x), repre- 
sented by thicker and thinner lines, respectively. They ap- 
pear inside the optical fiber when one of the two laser beams 
is turned off. We use r = 4/9. 

The limit r — * 1 describes the case of interactions be- 
tween laser beams of same amplitude but opposite circu- 
lar polarizations in an isotropic Kerr medium. This case 
is simpler, and can be described by a system of sym- 
metrically coupled nonlinear Schrodinger equations. For 
the pair of symmetrically coupled equations, if we set 
Fi(u,v) = Gi(u,v) = u 2 + v 2 we get the well known 
Manakov model p5[ . This model can also be used to 
represent the isotropic XY model |2(|, which obeys the 
Ginzburg-Landau equation dA/dt = d 2 A/dx 2 + A — 
\A\ 2 A, where A is complex. We use A = X + iY to show 
that for static configurations the above Ginzburg-Landau 
equation leads to the symmetrically coupled nonlinear 
Schrodinger equations. We use this as a guide toward 
the more general asymmetric model, and we consider the 
anisotropic XY model, which is described by 



OA 
~dt 



d 2 A 
dx 2 



A - \A\ 2 A + aA 



(14) 



where a is real. We use this model to investigate the case 
of two laser beams, in a medium such that F r (u,v) = 
u 2 + r(2r + l)v 2 and G r (u, v) = r 2 v 2 + r(2r + l)u 2 , that 
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obey F r (u, 0) 



and G r (0, v) = r v , as required to 



recuperate the case describing a single laser beam, when 
the other beam is turned off. The form of the interaction 
is dictated by the anisotropic XY model, and by the two 
laser beams that enter the fiber. 

The above considerations settle the problem. And now 
the key issue follows after recognizing that the coupled 
Schrodinger equations can be seen as equations of motion 
for static fields that appear in models of two real scalar 
fields in bidimensional space-time . In the case under 
examination the system of two fields is described by the 
potential 



V(u,v) = i (-1 + u 2 + rv 2 ) 2 + - (2ruvf 



This potential has the general form 28 3 



V{u,v) 



1 fdWy 1 fdW\ 

2 \~du~) + 2 \dv ) 



where W = W(u, v) is given by 



W(u, v) = —u + -u 3 + ruv 2 
o 

The pair of second order equations now become 



l a u 
2dx^ 
l d 2 v 
2dx^ 



-u + u 3 + r(2r + l)uv 2 
-rv + r 2 v 3 + r(2r + l)u 2 v 



(15) 
(16) 
(17) 

(18) 
(19) 



The main issue here is that the above system of nonlinear 
Schrodinger equations can be solved exactly. There are 
two trivial solutions, that appear when one turns off one 
of the two beams. They are the fundamental black soli- 
tons ui(x) = Ub(x), v i = and U2 = 0, V2{x) = Vb(x), 
which reproduce the former solutions (12) and (|l3|). Fur- 
thermore, we also have the non-trivial vector soliton, de- 
scribed by 



u 2{ x ) = i tanh[2r(a; — x)] 



i>2(x) = ±\/l/r — 2 sech[2r(a; — x)] 



(20) 
(21) 



This result is very interesting: it shows that two black 
solitons in interaction inside the silica fiber generate a 
pair of solitons, one of the black type and the other of 
the bright type. In Fig. 2 we depict this non-trivial vec- 
tor soliton. It illustrate that the presence of bright soli- 
tons is not a privilege of specific fibers, that have anoma- 
lous group-velocity dispersion. We may say that the sec- 
ond laser beam sees the normal group-velocity disper- 
sion of the medium as an anomalous one, and this ap- 
pear through the presence of the first laser beam, which 
is responsible for changing the group-velocity dispersion 
inside the fiber. This phenomenon is new, and follows 
from the vector soliton (^ti| ) and (fjl[). We notice that 
the vector soliton only exists for r G (0,1/2), and the 



limit r — > 1/2 leads us back to the scalar soliton Ub(x). 
This model has the property that the limit r — ► 1 is triv- 
ial. This limit leads to the symmetric situation, which 
can be investigated after rotating the the (u, v) plane to 
(u+,it_), with u± = u ± v. These transformations fac- 
torate the two-field system into two degenerate systems 
of a single field, that support no vector field solutions of 
the black and bright type. This shows that our model is 
asymmetric by construction, and the limit r — > 1 leads 
to a trivial symmetric model. 

The presence of the nontrivial black-bright vector soli- 
ton is easier to see if we integrate the coupled equations 
(|l8|) and dl^) once. This gives the first order equations 



1,2, 2 

— = =pl ± u ± rv 
dx 

dv 

— = ±2ruv 
dx 



(22) 
(23) 



The fact that the pair of second order equations ( |18[ ) and 
( |l9| ) can be integrated to give the above pair of first or- 
der equations (|22j) and ( |23|) is fundamental, since it is 
easier to integrate. In Field Theory this identifies the 
Bogomornyi-Prasad-Sommerfield bound [|TJ; the solu- 
tions are stable and minimize the energy of the system 

MM. 
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FIG. 2. The black (u|(x), thicker line) and bright {yl{x), 
thinner line) solitons that appear inside the optical fiber when 
both the laser beams are turned on. We use r = 4/9. 

We solve the first order equations using the trial or- 
bit method introduced in Ref. (3^]. We try the orbit 
v 2 + av 2 = 1, for a being some real constant. We 
differentiate this expression and we use the first order 
equations ( p2| ) and (^3|) to obtain u 2 + r(l + 2a)v 2 = 1. 
This is compatible with the previous choice if and only if 
a = r I (1 — 2r). Thus the orbit is given by 



1 



v 2 = l 



(24) 



It is an elliptical arc, that goes from [u = ±l,v = 0) to 
(u = =pl, v = 0) when x spans the entire real line. 

We use the orbit ( |24| ) to write the first order Eq. Q22 
in the form du/dx = ^2^(1 — u 2 ). It can be integrated 
easily. The result gives the solution (EOh. The other 
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solution ( pl| ) is obtained immediately from the orbit ( |24] ) 
obeyed by u and v. 

We have found the vector soliton as a solution of the 
pair of first order equations (0) and ( p3[) . The approach 
is new, and involves two steps: the first consists in trans- 
forming the pair of second order Schrodinger equations 
( |l8| ) and ( |i~9| ) into a pair of first order differential equa- 
tions; the second step deals with the trial orbit method 
introduced in Ref. p2| |. The trial orbit method was in- 
troduced to help solving second order equations of mo- 
tion for relativistic systems of coupled scalar fields. In 
the present context, the effectiveness of the trial orbit 
method is direct, since we can use the first order equa- 
tions themselves to check if the orbit we trial is good, 
that is, if the trial orbit does not contradict the first or- 
der equations that follow from the equations of motions. 
This check is immediate, and provides a direct trial for 
the eligibility of the orbit under investigation. This fact 
does not appear in the original proposal in Ref. p2f , since 
there it deals with second order differential equations. 

We see that when the potential is written in terms of 
W(u,v), in the form (|l6|), the two equations of motion 



are written as 



^2 = W uWuu + W V W VU 



dx 2 



w u w uv + w v w v 



(25) 
(26) 



where W u — dW/du, etc. We associate to these equa- 
tions the pair of first order equations 



dU TT/ dV TT/ 

-i-=W u , — = W V 

ax ax 



(27) 



We differentiate these first order equations to see that 
their solutions also solve the equations of motion. Thus, 
we can concentrate on the simpler task of solving first 
order equations to get solutions to the second order equa- 
tions. This is an important advantage, but it is restricted 
to work when the composite vector soliton presents non- 
trivial topological feature, as happens with dark solitons. 
This is not the case for bright solitons, that engender 
trivial topological behavior. Thus we may wander if it is 
possible to use two bright solitons to generate a compos- 
ite, bright-black vector soliton. Because the bright soli- 
tons are nontopological they cannot appear as solutions 
of first order equations, as recently shown in Ref. |}3| . 
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